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We identify what ideal correlated photon number states are to required to maximize the discrep- 
ancy between local realism and quantum mechanics when a quadrature homodyne phase measure- 
ment is used. Various Bell inequality tests are considered. 
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I. INTRODUCTION 

There has been recently active interest in tests of quan- 
tum mechanics [Q versus local realism in a high efficiency 
detection limit. Several authors including ourselves 
have considered detection schemes quadrature phase ho- 
modyne measurements. Such schemes use strong local 
oscillators and hence have very high detection efficiency 
||. This removes one of the current loopholes |^-^| and 
potentially allows a strong test of quantum mechanics 
JTof to be performed. 

The original idea of Gilchrist et. al. Q was to use a 
circle or pair coherent state |Tl] |f3|] produced by non- 
degenerate parametric oscillation with the pump mode 
adiabatically eliminated. Using highly efficient quadra- 
ture phase homodyne measurements, the Clauser Hornc 
strong Bell inequality jl4|-[l6| could be tested in an all 
optical regime. A small (approximately f .5%) but signif- 
icant theoretical violation was found for this extremely 
idea system. While the mean photon number for the 
system may be low (approximately 1.12), the use of ho- 
modyne measurements allow a macroscopic current to be 
detected. 

In this article, we take an unphysical but interesting 
approach and answer the following questions: 

• Given that your detection scheme is a quadrature 
phase homodyne measurement, what is the optimal 
input or correlated photon number state to maxi- 
mize the potential violation? 

• What is the optimal Bell inequality to test? 

To begin we will restrict our attention to correlated pho- 
ton number states of the form 



|*)=X) C „|n)|n) 



(1) 



n=0 



Two main sources of correlated photon number currently 
exist, each having it own particular form of c„. The most 
well known is simply the nondegenerate parametric am- 
plifier specified by an ideal Hamiltonian of the form Jl7[ 



where e is field amplitude of a nondepleting classical 
pump and \ is proportional to the susceptibility of the 
medium, a, b are the boson operators for the orthogonal 
signal and idler modes. After a time t, the state of the 
system is given by (0) with c„ specified by 



tanh™ [\et] 
cosh [%er] 



(3) 



In the quadrature phase amplitude basis this state has a 
positive Wigner function. Hence it can be described as 
a local hidden variable theory and thus cannot violate a 
Bell inequality. 

The other source of highly correlated photon number 
states exists in nondegenerate parametric oscillation. In 
the limit of very large parametric nonlinearity and high 
Q cavities, a state of the form [p3 
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can be generated. Here r is the size of the circle of the co- 
herent states and Jo is the zeroth order modified Bessel 
function. Equivalently this state can be written in the 
form of (GJ) with c„ given by 



(5) 
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H = -h X e (ab + a^tf) 



(2) 



This was the state considered by Gilchrist et. al. M 

Given the general form of known correlated number 
states ([!]), the next fundamental question that should be 
initially addressed is what we mean by the Bell inequality. 
A number of Bell inequalities exist, and the particular 
one depends used heavy on your application and exper- 
imental setup. The Bell inequalities to be considered in 
this article are the Clauser Home |lj| , the spin |lj] , and 
the information-theoretic Jl8| Bell inequality. A detailed 
derivation of the various inequalities will not be given, 
the reader is referred to references |l5| , |l^ , p^[ . Here we 
will consider only strong inequalities, that is inequalities 
where auxiliary assumptions (not based on local realism) 
are not required. In Fig (Q) we depict a very idealized 
setup for general Bell inequality experiment. 
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Probably the most well known inequality is the Clauser 
Home strong Bell inequality (l5| given by 



IBchl < 1 



(6) 



where 



B r 



Pn 



Pn 



■P ll (0,4f)+P 11 (0',<l/) 



Pi (9') + P 1 (0) 



(7) 



Here Pn is the probability that a "1" results occurs at 
each analyzer A, B given 9, <fi. Similarly Pi is the prob- 
ability that a "1" occurs at a detector while having no 
information about the second. For many of the actual 
experimental considerations an angle factorization occurs 
so that Pn (9, <f) depends only on 9 + 4>. Also Pi (9) and 
Pi (4>) are independent of 9, <f>. In this case B c h can be 
simplified to 



B c h 



3P n (VQ - Pn (3VQ 
2Pi 



(8) 



where ip = 9 + <j) = -9' - <p' = 9 + <p' and 3tp = 9' + <p. 

The second form of the Bell inequality (sometimes re- 
ferred to as the spin or original Bell inequality) is [jTjj 

B, = |£(M) -E{6',4>) 

+ E(0,<p')+E(e',<f>')\<2 (9) 



where the correlation function E I 



E(9, 



Pn 



Pool 
Poi 



is given by 



Poi(M) (io) 



Here as discussed above P n is probability that a "1" re- 
sults occurs at each analyzer A, B given 9, <f>. P o is prob- 
ability that a "0" results occurs at each analyzer A, B, 
while Pio (Poi) is probability that a "1" ("0") results oc- 
curs at the analyzer A and a "0" ("1") at B. With the 
angle factorization given above, the inequality (|^) can be 
rewritten as 



|3-E(V0 -E(3ip)\ < 2 



(11) 



Our final form of Bell inequality to be considered in 
this article was developed by Braunstein and Caves [Gl| . 
This classical information-theoretic Bell inequality has 
the form 



B 



info 



> o 



where 



B info = -H (9\<f>) +H{9\<t>') 

+ H{<t>'\6')+H{9'\<j>) 



Here H (0\<j>) is given by 



P(a,b) 
P(a) 



(12) 



(13) 



(14) 



with log (P (a, b)/P(a)) being the information gained at 
B given the result at A is known. The conditional in- 
formation is then given by H (9\(j>). The base of the log- 
arithm determines the units of the information (base 2 
for bits, base e for nats). For quantum computing pur- 
poses, this inequality should prove highly useful as it di- 
rectly deals with information content. Several other Bell 
inequality do exist such as the CHSH inequality ||, but 
these are not considered here due to there weaker nature. 
Auxiliary assumptions are necessary in there derivation 
which open up several loopholes |7|-p| . 



II. CORRELATED STATES 

From ([!]) we need to find the optimal c n which gives 
the largest Bell inequality violation. Before determin- 
ing the c„ we need to briefly focus our attention on the 
quadrature phase homodyne measurement. 

A quadrature phase-amplitude homodyne measure- 
ment X(0) at A can achieved by combining a signal field 
(say d) with a strong local oscillator field (say e) to form 
two new fields given by c± = [a ± eexp (iff)] /y/2. Here 
9 is a phase shift which allows the choice of particular 
observable to be measured, for instance choosing 9 as 
or 7r/2 allows the measurement of the conjugate phase 
variables X(Q) and X(ir/2) respectively. The homodyne 
measurement gives the photocurrent difference as 



Id = c\c+ - C_C- 



eX(9) 



(15) 



Performing a measurement on the quadrature phase am- 
plitude X(9) at A yields a result x\(9) which ranges in 
size and sign. Similarly a measurement on the quadra- 
ture phase amplitude X(d>) at B yields a result Xi (((>). 
For our state given by (|l|), the probability of obtaining 
the result xi(6),X2(<f) is simply 



p xlX2 (M)H<*i(0)i<*2Wi*>r 



(16) 



where 



{x(ip)\n) = 



v/2 



^e- x V*H n ( Xi ). (17) 



Here H n (xi) is the Hermite polynomial and tp is the phase 
of the local oscillator. Eqn (|l^) can be explicitly written 
as 



Px IK2 w = E 



3 — i(n— m)ip 
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(18) 
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where ip = 6 + (f>, that is our expression depends only on 
the sum of the individual local angles. 

The probability given by ( |l8| ) is for continuous vari- 
ables. The majority of the tests of quantum mechanics 
versus local realism require a binary result. Hence for 
a given quadrature measurement Xi we classify the re- 
sult as "1" if Xi > and the mutually exclusive "0" 
if Xi < 0. Here we have set the binning window about 
Xi = 0. Where this binning window is located is quite ar- 
bitrary, but the maximum violation occurs for the value 
we have selected. 

The probability of obtaining both particles in the "1" 
bin is 



poo pOO 

Pll W) = / / dx 1 dx 2 P xl x 
Jo Jo 



W) 



(19) 



while the probability of obtaining both particles in the 
"0" bin is 



Poo WO 







1) 



dx 1 dx 2 P xlX2 W) 



(20) 
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The other probabilities such Pio (tp) , Pqi (ip) can be cal- 
culated in a similar fashion. The probabilities formu- 
lated above are joint probabilities. Various of the strong 
Bell inequalities also require marginal probabilities of the 
form 



Pi W) 



OO pOO 
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(21) 



The above integrals can be easy evaluated using the 
results pi 



'H n (x)H m (x) 



7T2 



[J-"(n, m) — T{vn^ n)] 
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where F(n, m) is given by 



(for n 7^ to) 



(22) 



(23) 



with r being the Gamma function. Performing the inte- 
grals for (113) and (pfl) we find 



Pa WO = p oo wo 



2 n+m+1 irc n c* m 
n\m\(n — to) 2 



x [J-"(n, m) — !F(m, n)] cos [(n — m)ip] (24) 
Similarly Eqn ( |2l| ) simplifies to 

Pi = 1/2 (25) 



which is independent of the sum of the local oscillator 
angle ip. It is also simple to calculate the correlation 
function E(ip) 



z — ' n\m\(n — mY 



x cos \(n — m)ip] 



(26) 



Given the probabilities Pu, P 00 , ... it is also possible 
to calculate the conditional information H (0\4>) 



H(6\<t>) 



Pu log [2P n ] 
P w log [2P W ] 



P 00 log [2P 00 ] 
Poi log [2P i] 



(27) 



It is now possible to calculate the Clauser Home (g) 
and spin (^|) and information-theoretic (|l2|) Bell inequal- 
ities. Some insight into the problem can be achieved by 
a careful examination of the term 



■)n+m 



n\m\(n — m)' 



[J-(n, m) — J-(m, n)Y 



(28) 



which is present in all the joint probability distributions. 
This expression has several interesting features. First, 
as the difference between n and to becomes large, the 
smaller that the above expression contributes to any of 
the probability distributions. The main contribution for 
the expression comes from the case to = n ± 1. Second, 
when n — to is even, the above expression is zero. Fi- 
nally, as n is large the different between the n, m = n — 1 
and n + 1 , to = n elements for fixed large n vanishes and 
they reach an asymptotic limit which is smaller than the 
n = 1,to = case. If these higher order n terms domi- 
nate due to the choice of the c n in the probability formula, 
then the various Bell inequalities can not violated. This 
also has the implication that the mean photon number 
cannot be high if a violation is to occur and hence it is 
not a macroscopic test of quantum mechanics. 



III. A SIMPLE CASE 

To begin our investigations of the Bell inequalities, 
consider the case of we have only two photon pair states 
that is, 



|*)=co|0)|0)+c 1 |l)|l) 



(29) 



where for convenience we choose c n real. We also re- 
quire Cq + c\ — 1. The joint probability distributions are 
readily calculated and in fact 



Pn WO = Poo WO = 7 + — cos M 
PioW)=Po 1 W) = 7-— «»M 
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Calculating B c h and B s from (^]) and (g) we find 



(30) 

(31) 
(32) 
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B c h — — 
B M = 



C Ci 



2 7T 

4c ci 



{3 cos [tpo] — cos 



{3 cos [tpo] — cos [3ip]} 



Optimizing for the angle ip we find 

B c h 
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(33) 

(34) 
(35) 



(36) 

(37) 
(38) 



that is, \B c h\ < 1 and |£? s | < 2 for all cq. No viola- 
tion of the strong Clauser Home or spin Bell inequality 
is possible. 

For the information theoretic case we find 



Hty) = --log 







'1 + 2A" 




- A log 




1 - 2A 



where 



A 



2c ci 



cos [ip] 



(39) 



(40) 



The informational theoretic Bell inequality is given by 
Bmfo — 3H (ip) — H (Sip) > 0. A violation of this in- 
equality is possible if Bi n f a < 0. Unfortunately for all 
cqci and ip we have Bi„f > 0. 

No violation is possible for an y o f the Bell inequalities 
considered for the ideal state (£9[) when the detection 
scheme is based on homodyne quadrature phase mea- 
surements. If more correlated photon pairs are present 
can a violation be achieved? The obvious answer is yes, 
because of the recent work of Gilchrist et. al. The 
real question is how large this violation is? 



IV. NUMERICAL STUDIES 

Considering the expression (|l]) for the correlated pho- 
ton pairs, what are the optimal c„ coefficients to maxi- 
mize the violation. Because of the results indicated by 
Gilchrist et. al. and our previous discussion we antici- 
pate that the mean photon number per mode must be low 
to obtain a violation. Hence we will truncate the num- 
ber state basis at 10 (20) photon per mode. Performing 
a numerical optimization over all the c„, the optimal set 
is found to maximize the Clauser Horne and spin Bell 
inequality (Table ||). A plot of c„ versus n is depicted in 

Fig® ' 

It is interesting to now discuss some properties of these 
optimal c„. First the general shape of the c„ versus n 
curve shown in (^) is similar to that considered in the 
circle state by Gilchrist et. al. It is however not 
exactly the same (see (Table |). Given this optimal pa- 
rameter set, what is the maximum violation of the Bell 



inequalities we are considering. In Fig (||) we plot both 
the Clauser Horne and spin Bell inequalities versus ip. 

For the Clauser Horne Bell inequality the maximum 
violation corresponds to B c h = 1.019, while the maxi- 
mum violation for the spin Bell inequality corresponds 
to B s = 2.076. Interesting here is that the percentage 
violation of the spin inequality is approximately 3.8% 
compared with the 1.9% for the Clauser Horne case. This 
significantly increases the potential for an experiment to 
be performed provided such an experiment were not sig- 
nificantly more difficult. Also the results for the optimal 
c n set give a Clauser Horne Bell inequality violation that 
is approximately 20% greater than the circle state results 
of Gilchrist et. al {§. 

It is interesting to consider whether a greater viola- 
tion of the Bell inequality can be achieved with the state 
given by (§). To this end we show the effect of the varia- 
tion of both r and ip (sum of the local oscillator angles) 
for both the Clauser Horne and spin Bell inequalities in 
Fig ([|). As can be seen the spin Bell inequality can be 
violated far more significantly than the similar Clauser 
Horne case. In fact, as occurred previously the percent- 
age maximum violation in the spin inequality is twice 
that of the Clauser Horne result. 

In any of the analysis considered above we have not 
discussed errors, there sources and how they effect the 
potential violation. We will not present any significant 
details here in this article but refer the reader to || for 
such a decision. 

Our final Bell inequality to be considered is the Braun- 
stein and Caves |18| information-theoretic case. In 
Fig we plot Bi n j versus ip. No violation of the 
information-theoretic inequality is possible for any ip. 

A question here to be addressed is why two of the 
strong inequalities can be violated while this information- 
theoretic Bell inequality is far from being violated. In 
the binning process to give a binary result for a quadra- 
ture measurement, information must be discarded. The 
information-theoretic inequality is much more sensitive 
to this information loss than the Clauser Horne inequal- 
ity. Also why would we fundamentally expect all three 
inequalities to be violated. A violation of any of the in- 
equalities indicate a discrepancy between quantum me- 
chanics and local realism. 



V. CONCLUSION 

In this article we have place strict bounds on the opti- 
mal c„ coefficients for the state (Q) which maximizes the 
Clauser Horne and Spin Bell inequalities when a homo- 
dyne quadrature phase measurements is performed. The 
spin Bell inequality is violated by approximately 3.6% 
while the Clauser Horne inequality is violated by approx- 
imately 1.9%. The violation is small however due to the 
fact that we are discarding information in the binning 
process. In fact due to the information loss in the bin- 
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ning process the information theoretic Bell inequality is 
not violated in any regime. A larger violation cannot be 
obtained using homodyne measurements with the strong 
inequalities we have considered. 

While our optimal c„ coefficient give a slightly better 
violation than the pair coherent state, it is difficult to see 
how such a state could be generated. Closely examining 
the spin Bell inequality with the pair coherent state still 
indicates that a greater violation (approximately twice 
the size) is possible than for the other inequalities. This 
would make the test much more feasible provided the pair 
coherent state could be generated. In such a system the 
mean photon number is small, so this is not strictly a 
macroscopic test of quantum mechanics. It does however 
have a macroscopic nature due the strong local oscillator 
which means large photodetector currents are obtained. 

To conclude, quadrature phase homodyne measure- 
ment provide a mechanism for performing tests of the 
Bell inequality with highly efficient detection. This al- 
lows one of the loopholes in current experiments to be 
closed. However, due to the inherent information loss in 
the binning process, the violations are small but should 
be achievable. 
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FIG. 1. Schematic of a very generalized Bell experiment 



setup. After a source prepares two particles, these particles 
are directed out to the locations A and B. At each location 
there is an analyzer with an adjustable parameters 9, <f>. The 
particles are then detected resulting in a binary result "1" or 
"0" individually. These results can then be used to build up 
the statistics necessary to test the various Bell inequalities. 
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FIG. 2. Plot of c„ versus n 




FIG. 3. Plot of the Clauser Home (a) and Spin (b) Bell 
inequality versus tp. A violation occurs for the Clauser Home 
Bell inequality if B c h > 1. A violation of the spin Bell in- 
equality occurs for B s > 2. 
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FIG. 5. Plot of the information-theoretic Bell inequality 
versus ip. A violation is possible for B in f < 0. 




FIG. 4. Plot of the Clauser Home (a) and Spin (b) Bell 
inequality versus r and ip. A violation occurs for the CH Bell 
inequality if B c h > 1. A violation of the Spin Bell inequality 
occurs for B s > 2. 
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TABLE I. The optimal c n parameters to maximize the vi- 
olation of the Clauser Home and spin Bell inequalities. The 
c n values for the Gilchrist circle state are also given 
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Cn 


Eqn (|) with r ~ 1.12 





0.4990 


0.5495 


1 


0.6355 


0.6893 


2 


0.4760 


0.4323 


3 


0.3135 


0.1808 


4 


0.1465 


0.0567 


5 


0.0235 


0.0142 


6 


0.0075 


0.0029 


7 


0.0024 


0.0005 


B c h 


1.019 


1.016 


B a 


2.076 


2.064 



